General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



IW\SA 

Technical Memorandum 79678 


NIMBUS 6 DOPPLER PROCESSING 
USING THE FAIRBANKS 
CALIBRATION PLATFORM 


D. W. Koch and J. G. Lefler 

W79-16032 

Uncias 
13575 


(NASA-TM-79678) NIHBUS 6 DOPPLER PROCESSING 
USING THE FAIRBANKS CALIBRATION PLATFORM 

CSC! 22B 


(NASA) 18 p HC A02/MF A01 


G3/15 


OCTOBER 1978 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 




NIMBUS 6 DOPPLER PROCESSING USINQ THE 
FAIRBANKS CALIBRATION PLATFORM 

David W. Koch 
Joan G. Leflcr 

ABSTRACT 

The weighted least squares processing of Doppler shifted signals transmitted from Fairbanks 
to the NIMBUS 6 satellite is examined. The paper represents research conducted in support of the 
NASA satcilite aided Search and Rescue program. An estimated NIMBUS 6 ephemeris, accurate to 
1 .5 - 2.5 km and 0.5 - 2.5 m/s relative to a reference orhit, is obtained during the three day signal 
transmission period. This suggests updating the knowledge of the relay satellite ephemeris by one 
reference beacon is needed during the Search and Rescue demonstration. Per pass frequency offsets 
arc recoverable on the order of 600 - 610 Hz ±0.4 Hz. Residual statistics indicate an excelleni data 
fit with a near zero mean and an observed noise generally less than 2.0 Mz. 
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NIMBUS 6 DOPPLER PROCESSING USING THE 
FAIRBANKS CALIBRATION PLATFORM 


1. INTRODUCTION 

This paper examines the processing of some actual one-way Doppler shifted signals transmitted 
by the Fairbanks, Alaska calibration platform to the NIMBUS 6 meteorological satellite for three 
days during June 1977. Processing consisted of the weighted least squares estimation of the NIMBUS 
6 epocli state and frequency offsets for individual satellite passes over the platform. The accuracy of 
results is measured by a 6 day comparison of the estimated NIMBUS 6 ephemeris with a reference 
orbit of known accuracy and examination of residual statistics. 

Several considerations motivated the paper*. 

(1) to provide an example of the orbit determination accuracies obtainable for NASA appli- 
cations satellites using only uplink signals transmuted from fixed platforms. The use of 
reference beacons co-located with local user terminals (LUT) to update the relay satellite 
knowledge during Search and Rescue operations is under consideration (Ref. 1 ). 

(2) to observe the long term behavior of individual radial, along track, and cross track orbital 
position errors during both tracking and free propagation. Such behavior is of particular 
importance for the Search and Rescue demonstration onboard the Television Infra-Red 
Observation Satellite (TIROS-N) for which each orbital position component must satisfy 
a different accuracy level. 

(3) to extract from real data a measure of frequency offset. During Search and Rescue oper- 
ations the emergency b-^acon emits distress signals accurate to only several kHz. Conse- 
quently the offset must be extracted from the data tp remove its degrading influence 
upon beacon location estimates. 

The processing was performed using a weighted least squares parameter estimation program 
(Ref. 2-4). The program is coded in FORTRAN-H with double precision arithmetic and resides on 
permanent disk on an IBM 360/95 digital computer at Goddard Space Flight Center (GSFC), A 
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typical data reduction required 400Kbytesof core storage and tlirce minutes each of Central 
Processing Unit (CPU) and Inpiit/Output (I/O) time. 

2. DATA PROCESSING APPROACH 

A nominal NIMBUS 6 epliemeris from June 1 ~4, 1977 was generated from a typical epoch 
state selected from predicted vectors available at GSFC. The ephemeris was numerically integrated 
with variable stepsize using a power series solution tlirough 16tli order to the equations of motion. 

The Eartli’s gravitational field was modeled by the spherical harmonic coefficients from the 
Goddard Eartli Model-7 (GEM-7) (Ref. 5) truncated at degree and order 8. Constants for the 
Earth’s mass, radius, and flattening corresponding to GEM-7 were included. 

The calibraf’on platform at Fairbanks, Alaska nominally transmitted signals on a carrier fre- 
quency of 401.2 MHz ±5 kHz to the NIMBUS 6 meteorological satellite throughout June 1977 in 
conjunction with the Tropical Wind Energy Conversion and Reference Level Experiment (TWERLE) 
(Ref. 6). The platform is one of five worldwide employed to correct the NIMBUS 6 orbital along 
track error during daily operations. The transmitted signal experiences a frequency sliift upon re- 
ception by NIMBUS 6 because of the satellite’s motion relative to tlie platform. These Doppler 
shifted signals constituted the data type processed. Twelve passes of data scattered over the first 
three days of June 1977 were selected. These passes contained 140 data points collectively with an 
individual pass including between 8 and 15 points separated by one minute. Near horizon data 
points were ignored. 

A weighted least squares adjustment process was used to reduce these data. Software con- 
straints necessitated the conversion of the data from Doppler sliift to metric range rate by multi- 
plication with tile negative of tlie signal waveiengtii. Tlie signals were transmitted at 401.2 MHz on 
a waveiengtii of 0,75 meters. The eighteen adjusted parameters included the 6-dimensional NIMBUS 
6 epoch state and one frequency offset for each of the twelve passes. The a priori state chosen to 
initiate the iterative adjustment process was assumed poorly known. Individual components were 
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weighted with imcorrejn.tcd uncertainties of 10 km and 100 m/s. Tor cacii pass a zero value for the 
initial frequency offset with an uncertainty of 1 .3 kHz was postulated. 

Tile adjusted NIMBUS 6 epoch state was subsequently propagated forward for six days using 
modelling identical to that of tlic data processing. Tills propagated epliemeris was then differenced 
over this same interval from the most accurate available estimate of the true NIMBUS 6 epliemeris. 
This estimate was independently derived from Minitrack interferometer tracking data and has an 
advertised accuracy of 500 meters RMS. Tlius the errors in the estimate during both signal trans- 
mission and free propagation can be examined. 

3. RESULTS 

The primary parameters extracted from the Doppler shift data were the six components of the 
NIMBUS 6 epoch state. Figure 1 illustrates the successive corrections to the current position epoch 
state at each Iteration during the processing. Corrections to individual components (AX, AY, A7) 
and total position (Ar) are displayed. Corresponding corrections to the velocity vector arc found in 
Table 1 within the Appendix. 

The initial epoch state was adjusted by 7 km in total position after the first iteration. There- 
after the total corrections decrease rapidly to less than 30 meters after three iterations. Corrections 
to individual position components vary widely. After tlie third iteration the adjustment process 
stabilizes with only minor (i.e., less than 10 meters) corrections. Convergence occurs within five 
iterations to an estimated epoch state accurate to 412 meters and 38 cni/s RSS. Each velocity 
component is estimated to an equal degree of accuracy. 

Subsequent to adjustment, the estimated epoch state was propagated forward for six days 

under assumptions identical to those of the data processing. To obtain a measure of the long term 
accuracy of this propagated epliemeris, it was subtracted over this same time interval from the most 
accurate available estimate of the true NIMBUS 6 ephemeris. Figure 2 depicts the differences be- 
tween the two ephemerides and is an accurate measure of the error in the estimated NIMBUS 6 
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CORRECTIONS (m) 



Figure 1. Corrections to Current NIMBUS 6 Position Vector vs Iteration Number 

During Doppler Processing 
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Figure 2. Accuracy of NIMBUS 6 Epliemeris Estimate 



cphemcris. Tlie corresponding radial, along track, and cross track position errois are found in the 
Appendix. The 3 km level represents the TIROS-N orbital position accuracy requirement 
for the Search and Rescue demonstration. During the signal transmission period from June 1 - 3, 
both position and velocity errors arc periodic. The former oscillates between 1 .5 and 2,5 km and 
the latter between 0.5 and 2.5 m/s. Their complementary beiiavior typifies the conservation of 
angular momentum of the NIMBUS 6 orbit. Outside the transmission period both errors exhibit 
continued oscillation superimposed upon a strong secular growth. Note the required TIROS-N ac- 
curacy requirement is easily achieved but only during periods of signal transmission. Tlius during 
tlie Searcli and Rescue demonstration continuous updating by at least one reference beacon is 
necessary. 

The transmitter at the Fairbanks, Alaska calibration platform has a nominal frequency of 
401 .2 MHz ±5 kHz. In reality, the true transmit frequency deviates continuously from this nominal. 
Deviations occur even over very sliort time periods sucli as tlie 15 minutes required for a typical 
satellite pass over a transmitter. They arc caused by crystal oscillator instability and normal degra- 
dation of iiardware components. The deviation between the true and nominal frequency is termed 
frequency offset. During Search and Rescue operations, an emergency beacon will transmit signals 
accurate to only several kHz. Thus the offset must be adjusted to remove its degrading effect upon 
the more important beacon location estimates. 

Table I lists estimates of tlie frequency offsets witli their accuracies for each of tlie 12 passes 
of NIMBUS 6 over the Fairbanks platform during June 1-3, 1977, A remarkable uniformity in 
.magnitudes of estimates is evident. All offsets are approximately recoverable between 600 to 
610 Hz accurate to 0.4 Hz. Estimates are independent of the number of data points processed. 

The accuracy of the estimated NIMBUS 6 ephemeris was determined previously by comparison 
with a reference orbit (Fig. 2). Another technique commonly employed to assess the accuracy of 
estimates is examination of data residuals. A parameter estimation program uses geometry to gen- 
erate data. These internally generated data, i.e. the “computed” measurements, are then 
compared point by point to the data being processed, 
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Table 1 

Estimates of Frequency Offsets and their Accuracy for 
NIMBUS 6 Passes over Fairbanks, Alaska Calibration Platform June 1 - 3, 1977 


Pass 

No. 


No. of 
Data Points 

Estimates 



1 

9 

9 

607 

0.4 

2 

10 

8 

611 

0.5 

3 

12 

12 

608 

0.4 

4 

13 

13 

*601 

0.4 

5 

11 

12 

608 

0.4 

6 

10 

9 

611 

0.4 

7 

12 

11 

605 

0.4 

8 

14 

13 

605 

0.4 

9 

11 

12 

601 

0.4 

10 

13 

14 

604 

0.4 

11 

14 

15 

596 

0.4 

12 

9 

10 

595 

0.4 


i.e. the “observed” measurements. The differences between the two are labelled “residuals” and 
their magnitudes reflect the quality of tiie data fit. The smaller tlie residuals implies tlie better the 
fit. Table 2 lists the mean value ifi) and standard deviation (a) of the Doppler shift residuals for the 
12 NIMBUS 6 passes over the Fairbanks platform. The near zero mean value indicates an excellent 
least squares fit to the data. The magnitude of this statistic is relatively independent of the number 
of data points. The standard deviation is a measure of the noise on the actual data. Table 2 reveals 
their magnitudes are generally less than 2 Hz, a value consistent with TWERLE frequency resolution 
of ±1 Hz, 
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Table 2 


Mean Value (/i) and Standard Deviation (a) of Doppler Shift Residuals for 
NIMBUS 6 Passes over Fairbanks, Alaska Calibration Platform June 1 - 3, 1977 


Pass 

No. 


No. of 
Points 

fX*‘ 

(10-'* Hz) 

o** 

(Hz) 

1 

9 

9 

0,68 

1.60 

2 

10 

8 

0.77 

1.84 

3 

12 

12 

0.51 

0.77 

4 

13 

13 

0.47 

1.35 

. 5 

11 

12 

0.52 

0.72 

6 

10 

9 

0,70 

1.29* 

7 

12 

1! 

0.57 

1.03 

8 

14 


0.48 

2.00 

9 

11 

12 

0.54 

1.33 

10 

13 

14 

0.48 

1.20 

11 

14 

15 

0.46 

1.27 

12 

9 

10 

0,67 

0.97 


"'estimate of quality of fit 
"""estimate of data noise 

4. CONCLUSIONS 


This paper has presented results of wciglited least squares processing of actual NIMBUS 6 Dop- 
pler shift data resulting from transmission by the Fairbanks, Alaska calibration platform on June 
1 - 3, 1977. The paper was prompted primarily by considerations arising from satellite aided Search 
and Rescue operations. Results indicate that a NIMBUS 6 ephemeris accurate to 1 .5 - 2,5 km in po- 
sition and 0.5 - 2.5 m/s in velocity is possible during periods of signal transmission. Outside this 
period both errors grow without bound. This implies that updating by one reference beacon is 
needed to meet orbital accuracy requirements during the Search and Rescue demonstration onboard 
TIROS-N. The frequency of an emergency transmitter will usually only be known to several kHz. 
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Hence position location software must estimntc frequency offsets in addition to tlie emergency 

transmitter location. This paper ha's ?:nown that offsets are recoverable on the order of 600 -CIO Hz 

to an accuracy of 0.4 Hz. Finally, statistics of the measurement residuals indicate an excellent least 

squares fit to the Doppler data and an observed noise consistent with expected values. 
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APPENDIX 


Tabic A-1 

CorrcctloDS to Current NIMBUS 6 Velocity Vector vs 
Iteration Number During Doppler Processing 


1 

Iteration 

Number 

(cm/s) 

AY 

(cm/s) 

bJL 

(cm/s) 

AV 

(cm/s) 

1 

-356 

-423 

-355 

657 

2 

+4 

-1 

+18 

18 

3 

-1 

-2 

-1 

2 

4,5 

0 

0 

0 

0 
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Figure A-2. Accuracy of NIMBUS 6 Along Track Position Estimate 
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Figure A-3. Accuracy of NIMBUS 6 Cross Track Po-'tion Estimate 



